Copy number variation and transposable elements feature in recent, ongoing adaptation at the Cyp6g1 locus

The increased transcription of the Cyp6g1 gene of Drosophila melanogaster, and consequent resistance to insecticides such as DDT, is a widely cited example of adaptation mediated by cis-regulatory change. A fragment of an Accord transposable element inserted upstream of the Cyp6g1 gene is causally associated with resistance and has spread to high frequencies in populations around the world since the 1940s. Here we report the existence of a natural allelic series at this locus of D. melanogaster, involving copy number variation of Cyp6g1, and two additional transposable element insertions (a P and an HMS-Beagle). We provide evidence that this genetic variation underpins phenotypic variation, as the more derived the allele, the greater the level of DDT resistance. Tracking the spatial and temporal patterns of allele frequency changes indicates that the multiple steps of the allelic series are adaptive. Further, a DDT association study shows that the most resistant allele, Cyp6g1-[BP], is greatly enriched in the top 5% of the phenotypic distribution and accounts for ~16% of the underlying phenotypic variation in resistance to DDT. In contrast, copy number variation for another candidate resistance gene, Cyp12d1, is not associated with resistance. Thus the Cyp6g1 locus is a major contributor to DDT resistance in field populations, and evolution at this locus features multiple adaptive steps occurring in rapid succession.

Analysis of genome-wide association study data using the protein knowledge base

Article 98

Bone marrow chimeric mice reveal a role for CX3CR1 in maintenance of the monocyte-derived cell population in the olfactory neuroepithelium

Macrophages in the olfactory neuroepithelium are thought to play major roles in tissue homeostasis and repair. However, little information is available at present about possible heterogeneity of these monocyte-derived cells, their turnover rates, and the role of chemokine receptors in this process. To start addressing these issues, this study used Cx3cr1gfp mice, in which the gene sequence for eGFP was knocked into the CX3CR1 gene locus in the mutant allele. Using neuroepithelial whole-mounts from Cx3cr1gfp/+ mice, we show that eGFP+ cells of monocytic origin are distributed in a loose network throughout this tissue and can be subdivided further into two immunophenotypically distinct subsets based on MHC-II glycoprotein expression. BM chimeric mice were created using Cx3cr1gfp/+ donors to investigate turnover of macrophages (and other monocyte-derived cells) in the olfactory neuroepithelium. Our data indicate that the monocyte-derived cell population in the olfactory neuroepithelium is actively replenished by circulating monocytes and under the experimental conditions, completely turned over within 6 months. Transplantation of Cx3cr1gfp/gfp (i.e., CX3CR1-deficient) BM partially impaired the replenishment process and resulted in an overall decline of the total monocyte-derived cell number in the olfactory epithelium. Interestingly, replenishment of the CD68lowMHC-II+ subset appeared minimally affected by CX3CR1 deficiency. Taken together, the established baseline data about heterogeneity of monocyte-derived cells, their replenishment rates, and the role of CX3CR1 provide a solid basis to further examine the importance of different monocyte subsets for neuroregeneration at this unique frontier with the external environment.

GAD2 on chromosome 10p12 is a candidate gene for human obesity

The gene GAD2 encoding the glutamic acid decarboxylase enzyme (GAD65) is a positional candidate gene for obesity on Chromosome 10p11–12, a susceptibility locus for morbid obesity in four independent ethnic populations. GAD65 catalyzes the formation of γ-aminobutyric acid (GABA), which interacts with neuropeptide Y in the paraventricular nucleus to contribute to stimulate food intake. A case-control study (575 morbidly obese and 646 control subjects) analyzing GAD2 variants identified both a protective haplotype, including the most frequent alleles of single nucleotide polymorphisms (SNPs) +61450 C>A and +83897 T>A (OR = 0.81, 95% CI [0.681–0.972], p = 0.0049) and an at-risk SNP (−243 A>G) for morbid obesity (OR = 1.3, 95% CI [1.053–1.585], p = 0.014). Furthermore, familial-based analyses confirmed the association with the obesity of SNP +61450 C>A and +83897 T>A haplotype (χ2 = 7.637, p = 0.02). In the murine insulinoma cell line βTC3, the G at-risk allele of SNP −243 A>G increased six times GAD2 promoter activity (p < 0.0001) and induced a 6-fold higher affinity for nuclear extracts. The −243 A>G SNP was associated with higher hunger scores (p = 0.007) and disinhibition scores (p = 0.028), as assessed by the Stunkard Three-Factor Eating Questionnaire. As GAD2 is highly expressed in pancreatic β cells, we analyzed GAD65 antibody level as a marker of β-cell activity and of insulin secretion. In the control group, −243 A>G, +61450 C>A, and +83897 T>A SNPs were associated with lower GAD65 autoantibody levels (p values of 0.003, 0.047, and 0.006, respectively). SNP +83897 T>A was associated with lower fasting insulin and insulin secretion, as assessed by the HOMA-B% homeostasis model of β-cell function (p = 0.009 and 0.01, respectively). These data support the hypothesis of the orexigenic effect of GABA in humans and of a contribution of genes involved in GABA metabolism in the modulation of food intake and in the development of morbid obesity.

Adamts5, the gene encoding a proteoglycan-degrading metalloprotease, is expressed by specific cell lineages during mouse embryonic development and in adult tissues.

The secreted metalloprotease ADAMTS5 is implicated in destruction of the cartilage proteoglycan aggrecan in arthritis, but its physiological functions are unknown. Its expression profile during embryogenesis and in adult tissues is therefore of considerable interest. β-Galactosidase (β-gal) histochemistry, enabled by a LacZ cassette inserted in the Adamts5 locus, and validated by in situ hybridization with an Adamts5 cRNA probe and ADAMTS5 immunohistochemistry, was used to profile Adamts5 expression during mouse embryogenesis and in adult mouse tissues. Embryonic expression was scarce prior to 11.5 days of gestation (E11.5) and noted only in the floor plate of the developing brain at E9.5. After E11.5 there was continued expression in brain, especially in the choroid plexus, peripheral nerves, dorsal root ganglia, cranial nerve ganglia, spinal and cranial nerves, and neural plexuses of the gut. In addition to nerves, developing limbs have Adamts5 expression in skeletal muscle (from E13.5), tendons (from E16.5), and inter-digital mesenchyme of the developing autopod (E13.5–15.5). In adult tissues, there is constitutive Adamts5 expression in arterial smooth muscle cells, mesothelium lining the peritoneal, pericardial and pleural cavities, smooth muscle cells in bronchi and pancreatic ducts, glomerular mesangial cells in the kidney, dorsal root ganglia, and in Schwann cells of the peripheral and autonomic nervous system. Expression of Adamts5 during neuromuscular development and in smooth muscle cells coincides with the broadly distributed proteoglycan versican, an ADAMTS5 substrate. These observations suggest the major contexts in which developmental and physiological roles could be sought for this protease.

Diagnosis of {alpha}+-thalassemias by determining the ratio of the two {alpha}-globin gene copies by oligonucleotide hybridization and melting curve analysis

Most cases of {alpha}-thalassemia result from large deletions at the {alpha}-globin locus (1). The {alpha}-globin gene cluster contains a tandem array of 2 nearly identical {alpha}-globin genes (HBA; Fig. 1A ) (2). The {alpha}0-thalassemias are characterized by deletions that inactivate both {alpha}-globin genes of a given chromosome, whereas in {alpha}+-thalassemias, one gene remains functional. The most widespread {alpha}+-thalassemias are those designated –{alpha}3.7 and –{alpha}4.2, according to the lengths of the deleted fragments (3).

Characterisation of monotreme caseins reveals lineage-specific expansion of an ancestral casein locus in mammals

Using a milk-cell cDNA sequencing approach we characterised milk-protein sequences from two monotreme species, platypus (Ornithorhynchus anatinus) and echidna (Tachyglossus aculeatus) and found a full set of caseins and casein variants. The genomic organisation of the platypus casein locus is compared with other mammalian genomes, including the marsupial opossum and several eutherians. Physical linkage of casein genes has been seen in the casein loci of all mammalian genomes examined and we confirm that this is also observed in platypus. However, we show that a recent duplication of β-casein occurred in the monotreme lineage, as opposed to more ancient duplications of α-casein in the eutherian lineage, while marsupials possess only single copies of α- and β-caseins. Despite this variability, the close proximity of the main α- and β-casein genes in an inverted tail-tail orientation and the relative orientation of the more distant kappa-casein genes are similar in all mammalian genome sequences so far available. Overall, the conservation of the genomic organisation of the caseins indicates the early, pre-monotreme development of the fundamental role of caseins during lactation. In contrast, the lineage-specific gene duplications that have occurred within the casein locus of monotremes and eutherians but not marsupials, which may have lost part of the ancestral casein locus, emphasises the independent selection on milk provision strategies to the young, most likely linked to different developmental strategies. The monotremes therefore provide insight into the ancestral drivers for lactation and how these have adapted in different lineages.

Variation in the gene coding for the M5 muscarinic receptor (CHRM5) influences cigarette dose but is not associated with dependence to drugs of addiction : evidence from a prospective population based cohort study of young adults

Background: The mesolimbic structures of the brain are important in the anticipation and perception of reward. Moreover, many drugs of addiction elicit their response in these structures. The M5 muscarinic receptor (M5R) is expressed in dopamine-containing neurones of the substantia nigra pars compacta and ventral tegmental area, and regulates the release of mesolimbic dopamine. Mice lacking M5R show a substantial reduction in both reward and withdrawal responses to morphine and cocaine. The CHRM5, the gene that codes for the M5R, is a strong biological candidate for a role in human addiction. We screened the coding and core promoter sequences of CHRM5 using denaturing high performance liquid chromatography to identify common polymorphisms. Additional polymorphisms within the coding and core promoter regions that were identified through dbSNP were validated in the test population. We investigated whether these polymorphisms influence substance dependence and dose in a cohort of 1947 young Australians.

Results: Analysis was performed on 815 participants of European ancestry who were interviewed at wave 8 of the cohort study and provided DNA. We observed a 26.8% increase in cigarette consumption in carriers of the rs7162140 T-allele, equating to 20.1 cigarettes per week (p=0.01). Carriers of the rs7162140 T-allele were also found to have nearly a 3-fold increased risk of developing cannabis dependence (OR=2.9 (95%CI 1.1-7.4); p=0.03).

Conclusion: Our data suggest that variation within the CHRM5 locus may play an important role in tobacco and cannabis but not alcohol addiction in European ancestry populations. This is the first study to show an association between CHRM5 and substance use in humans. These data support the further investigation of this gene as a risk factor in substance use and dependence.

Candidate disease gene prediction using Gentrepid: application to a genome-wide association study on coronary artery disease

Current single-locus-based analyses and candidate disease gene prediction methodologies used in genome-wide association studies (GWAS) do not capitalize on the wealth of the underlying genetic data, nor functional data available from molecular biology. Here, we analyzed GWAS data from the Wellcome Trust Case Control Consortium (WTCCC) on coronary artery disease (CAD). Gentrepid uses a multiple-locus-based approach, drawing on protein pathway- or domain-based data to make predictions. Known disease genes may be used as additional information (seeded method) or predictions can be based entirely on GWAS single nucleotide polymorphisms (SNPs) (ab initio method). We looked in detail at specific predictions made by Gentrepid for CAD and compared these with known genetic data and the scientific literature. Gentrepid was able to extract known disease genes from the candidate search space and predict plausible novel disease genes from both known and novel WTCCC-implicated loci. The disease gene candidates are consistent with known biological information. The results demonstrate that this computational approach is feasible and a valuable discovery tool for geneticists.

Tumor progression locus 2 (Tpl2) deficiency does not protect against obesity-induced metabolic disease

Obesity is associated with a state of chronic low grade inflammation that plays an important role in the development of insulin resistance. Tumor progression locus 2 (Tpl2) is a serine/threonine mitogen activated protein kinase kinase kinase (MAP3K) involved in regulating responses to specific inflammatory stimuli. Here we have used mice lacking Tpl2 to examine its role in obesity-associated insulin resistance. Wild type (wt) and tpl22/2 mice accumulated comparable amounts of fat and lean mass when fed either a standard chow diet or two different high fat (HF) diets containing either 42% or 59% of energy content derived from fat. No differences in glucose tolerance were observed between wt and tpl22/2 mice on any of these diets. Insulin tolerance was similar on both standard chow and 42% HF diets, but was slightly impaired in tpl22/2 mice fed the 59% HFD. While gene expression markers of macrophage recruitment and inflammation were increased in the white adipose tissue of HF fed mice compared with standard chow fed mice, no differences were observed between wt and tpl2 mice. Finally, a HF diet did not increase Tpl2 expression nor did it activate Extracellular Signal-Regulated Kinase 1/2 (ERK1/2), the MAPK downstream of Tpl2. These findings argue that Tpl2 does not play a non-redundant role in obesity associated metabolic dysfunction.

Is there evidence of selection in the dopamine receptor D4 gene in Australian invasive starling populations?

 Although population genetic theory is largely based on the premise that loci under study are selectively neutral, it has been acknowledged that the study of DNA sequence data under the influence of selection can be useful. In some circumstances, these loci show increased population differentiation and gene diversity. Highly polymorphic loci may be especially useful when studying populations having low levels of diversity overall, such as is often the case with threatened or newly established invasive populations. Using common starlings Sturnus vulgaris sampled from invasive Australian populations, we investigated sequence data of the dopamine receptor D4 gene (DRD4), a locus suspected to be under selection for novelty-seeking behaviour in a range of taxa including humans and passerine birds. We hypothesised that such behaviour may be advantageous when species encounter novel environments, such as during invasion. In addition to analyses to detect the presence of selection, we also estimated population differentiation and gene diversity using DRD4 data and compared these estimates to those from microsatellite and mitochondrial DNA sequence data, using the same individuals. We found little evidence for selection on DRD4 in starlings. However, we did find elevated levels of within-population gene diversity when compared to microsatellites and mitochondrial DNA sequence, as well as a greater degree of population differentiation. We suggest that sequence data from putatively nonneutral loci are a useful addition to studies of invasive populations, where low genetic variability is expected.

Variation at the DRD4 locus is associated with wariness and local site selection in urban black swans

 BACKGROUND: Interactions between wildlife and humans are increasing. Urban animals are often less wary of humans than their non-urban counterparts, which could be explained by habituation, adaptation or local site selection. Under local site selection, individuals that are less tolerant of humans are less likely to settle in urban areas. However, there is little evidence for such temperament-based site selection, and even less is known about its underlying genetic basis. We tested whether site selection in urban and non-urban habitats by black swans (Cygnus atratus) was associated with polymorphisms in two genes linked to fear in animals, the dopamine receptor D4 (DRD4) and serotonin transporter (SERT) genes.

RESULTS: Wariness in swans was highly repeatable between disturbance events (repeatability = 0.61) and non-urban swans initiated escape from humans earlier than urban swans. We found no inter-individual variation in the SERT gene, but identified five DRD4 genotypes and an association between DRD4 genotype and wariness. Individuals possessing the most common DRD4 genotype were less wary than individuals possessing rarer genotypes. As predicted by the local site selection hypothesis, genotypes associated with wary behaviour were over three times more frequent at the non-urban site. This resulted in moderate population differentiation at DRD4 (FST = 0.080), despite the sites being separated by only 30 km, a short distance for this highly-mobile species. Low population differentiation at neutrally-selected microsatellite loci and the likely occasional migration of swans between the populations reduces the likelihood of local site adaptations.

CONCLUSION: Our results suggest that wariness in swans is partly genetically-determined and that wary swans settle in less-disturbed areas. More generally, our findings suggest that site-specific management strategies may be necessary that consider the temperament of local animals.

Tripartite motif-containing 55 identified as functional candidate for spontaneous cardiac hypertrophy in the rat locus cardiac mass 22

BACKGROUND: Left ventricular (LV) hypertrophy is a risk factor for cardiovascular death, but the genetic factors determining LV size and predisposition to hypertrophy are not well understood. We have previously linked the quantitative trait locus cardiac mass 22 (Cm22) on chromosome 2 with cardiac hypertrophy independent of blood pressure in the spontaneously hypertensive rat. From an original cross of spontaneously hypertensive rat with F344 rats, we derived a normotensive polygenic model of spontaneous cardiac hypertrophy, the hypertrophic heart rat (HHR) and its control strain, the normal heart rat (NHR).

METHODS AND RESULTS: To identify the genes and molecular mechanisms underlying spontaneous LV hypertrophy we sequenced the HHR genome with special focus on quantitative trait locus Cm22. For correlative analyses of function, we measured global RNA transcripts in LV of neonatal HHR and NHR and 198 neonatal rats of an HHR × NHR F2 crossbred population. Only one gene within locus Cm22 was differentially expressed in the parental generation: tripartite motif-containing 55 (Trim55), with mRNA downregulation in HHR (P < 0.05) and reduced protein expression. Trim55 mRNA levels were negatively correlated with LV mass in the F2 cross (r = -0.16, P = 0.025). In exon nine of Trim55 in HHR, we found one missense mutation that functionally alters protein structure. This mutation was strongly associated with Trim55 mRNA expression in F2 rats (F = 10.35, P < 0.0001). Similarly, in humans, we found reduced Trim55 expression in hearts of subjects with idiopathic dilated cardiomyopathy.

CONCLUSION: Our study suggests that the Trim55 gene, located in Cm22, is a novel candidate gene for polygenic LV hypertrophy independent of blood pressure.